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Abstract 

We probe electro-mechanical properties of InAs nanowire (diameter ~100 nm) resonators where 
the suspended nanowire (NW) is also the active channel of a field effect transistor (FET). We 
observe and explain the non-monotonic dispersion of the resonant frequency with DC gate voltage 
iy^^)- The effect of electronic screening on the properties of the resonator can be seen in the 
amplitude. We observe the mixing of mechanical modes with ■ We also experimentally 
probe and quantitatively explain the hysteretic non-linear properties, as a function of V^'^ , of the 
resonator using the Duffing equation. 
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Nano electro mechanical systems (NEMS) are being used extensively to study small 



ricated using carbon nanotubes 







12h16[. graphene 



displacements [2| , mass sensing |3|-|6| , spin-torque effect j3] , charge sensing j8( , Casimir force 
9| and potential quantum mechanical devices 



h 



.11 . A variety of NEMS devices, fab- 
ijisl' nanowires (NWs) of silicon 20 1 
and by micromachining bulk silicon 2l|], have been used to probe the underlying physics 
at nano scale. In this work we study the electromechanical properties of doubly clamped 
suspended n-type InAs NWs. In our suspended NW FET, the gate electrode serves three 
purposes: first, to modify the tension in the NW, second, to actuate the mechanical motion 
of the resonator and third, enabling us to systematically study the coupling of mechanical 
properties to the tunable electron density. As we will show, a tunable electron density leads 
to a variable screening length of the order of the nanowire's cross sectional dimensions. Thus 
the electro-mechanical properties enter into a mesoscale regime. Such a variable electron 
density is not accessible in carbon nanotubes; the screening length cannot be tuned contin- 
uously - relative to the diameter of the carbon nanotube - as easily. However, the physics 



of charge screening in nanoscale capacitors 22|] and ferroelectric devices 23j is intimately 
connected to that in our NEMS devices. Taken together, these observations suggest that 
a geometrical interpretation of capacitance is inadequate at the nano scale. Additionally, 
the gate allows us to tune the resonant frequency non-monotonically due to the competition 
between the electrostatic force and the mechanical stiffness (~ 1 N/m) of the nanowire, 
a feature expected (but heretofore not studied in detail) for all electrostatically actuated 
NEMS. In this article, we demonstrate that the InAs semiconducting nanowire system man- 
ifests this transition from softening to hardening as the gate voltage is varied. In addition, 
mixing of the natural mechanical modes as a function of V^^*" can be understood in terms 
of the structural asymmetries in the resonator. In the non-linear regime we study in de- 
tail hysteretic behavior as a function of V^'" (unlike the commonly studied response as a 
function of drive frequency) and we show that this can be understood by using the Duffing 
equation incorporating the effect of gate voltage. The observed hysteretic response with 
V^'" is an alternate knob for tuning the nonlinear response of our NEMS devices and can 
be used for charge detection 



2J]. Our work provides further understanding of the unique 



characteristics of NEMS devices operating at room temperature. The observed behavior can 
provide information on the nanomechanics of other systems whose electron density, stiffness 
or screening length cannot be so readily tuned. 
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The InAs NWs used for this work were grown using metal organic vapor phase epitaxy 



(MOVPE) |25|,|26|. The NWs are oriented in the < 111 > direction and are 80 to 120 nm in 
diameter with a length of several micrometers. The substrate used for making the devices is a 
degenerately doped silicon wafer with 300 nm thick Si02. We have fabricated suspended InAs 
NW devices by sandwiching the NWs between two layers of electron-beam resist and then 
using electron beam lithography to define the electrodes and suspend them by depositing 
~150 nm Cr and ~250 nm Au after development followed by in situ plasma cleaning 26 |. 
The ohmic contacts also serve as mechanical supports for the NW suspended ~200 nm 
above the surface of Si02. FigJT^ shows the SEM image of a resonator device and scheme 
for actuating and detecting the motion of the resonator. All the measurements were done 
at 300 K and pressure less than 1 mBar. 

To actuate and detect the resonance we use the device as a heterodyne mixer 
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29|. We use electrostatic interaction between the wire and gate to actuate the motion in a 
plane perpendicular to the substrate. We apply a radio frequency (RF) signal Vg{u) and a 
DC voltage V^'^ at the gate terminal using bias-tee. Another RF signal Vsd{^ + Aw) is 
applied to the source electrode (FigH^). The RF signal applied at the gate Vg{u) modulates 
the gap between NW and substrate at angular frequency u, and V^^*" alters the overall 
tension in the NW. The amplitude of the current through the NW at the difference frequency 
(Au), also called the mixing current Imix{Au), can be written as jl2l | 

(Au) = Ioscn{Au) + hack{Auj) 

where G is the conductance of the NW, q is the charge induced by the gate, Cg 
is the gate capacitance, z{uj) is the amplitude of oscillation at the driving frequency uj 
and z axis is perpendicular to the substrate. The term Ij,ack{Au) = ^^CgVg{u)VsD is 
the background mixing current which is independent of the oscillation of the NW and 
Ioscn{Au) = ^^^^z{u)V^'-'VsD depends on the amplitude of oscillation. Measuring 
Ifnix{Au}) using a lock- in allows us to monitor the resonance of the NW as the frequency 
is swept. FigHb shows ImixiAu) as a function of u for V^^"" = ± 36.4 V. The sharp fea- 
ture corresponds to the mechanical resonance of the NW. We address the asymmetry of 
the mixing current signal for iV^^*" later. FigHJ:; shows the plot of conductance (G) as a 
function of . The variation of G with Vg^'" is very critical for this scheme of heterodjTie 



mixing to work as it controls the overall amplitude Imix{^^)- The hysteresis observed in 
the measurement of conductance is typical for our suspended devices and is associated with 



charge trap states with dipolar nature on surface of the nanowire 30|, l31 |. 



We can connect the resonant frequency of the fundamental mode, /o, of a doubly clamped 
beam at zero V/^*", to the material properties of the beam as /o = Cqt^K where Vg — 



is the velocity of sound; r is the radius of the beam, / is the length of the beam, E is the 
Young's modulus, p is the density of the material and Cq = 1.78. Fig{T]i shows a plot for Vs 
that does not vary much from the bulk value 32] (dashed line) for eight different devices. 
The scatter around the Vg calculated using bulk values could be due to the relative volume 
fraction contribution of the amorphous layer around the NWs; this needs further detailed 
study. 

Figj2^ and 2b show the colorscale plot of Imix{^^) as a function of V^'~^ and a; on a 
logscale spanning more than 3 decades. Data is taken by sweeping V^'~^ for each value of 
frequency. In the data from both the devices we see a symmetric evolution of the resonant 
frequency as a function of . The parabolic behavior is expected as the attractive force 
exerted by the gate on the wire is given by Foe = \iy^^Y^C!g. An increase in V^'~' 



enhances the tension in the NW 



33|. We now discuss the particular W-shaped dispersion of 



modes as a function of V^^ . As \V^^\ is increased, initially the mode disperses negatively 
and after a certain threshold voltage Vg^ it disperses with a positive slope. Although purely 
negatively and positively dispersing modes have been studied in detail before by Kozinsky et 



al. 



3J], we observe the crossover regime where these interactions compete. This particular 



dispersion can be understood using a toy-model in which a wire is suspended from a spring 
of force constant Kj above the substrate. The wire and the substrate make up the two 
electrodes of a capacitor. There are two consequences of increasing first, it changes 

the equihbrium position by moving the NW closer to the substrate and second, it makes 
the local potential asymmetric and less steep. The effective force constant, Kgjj, is reduced, 
resulting in negative dispersion of the mode for iVg^*"! < |V^*^|. If the intrinsic force constant 
Kj had been independent of V^^*" , the modes would always disperse negatively when motion 
occurs perpendicular to the gate plane. However, in general Kj can be written as + 
aiy^cy ^ Piy^^Y + H.O.iVg^^)); where k, a, and /3 are constants. With V^^^ ^ 
then, to a first approximation. Kg// = Kj — ^{yf^Y^^^^j^- We find that the W-shape of the 
dispersion curve can be explained only if one considers the case where /3 > 0. The result of 
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such calculations (Figj2]:) quantitatively explains the experimental observations. Vg^, the 
value at which the crossover from negative to positive dispersion occurs is a function of the 
dimensions of the NW and the capacitive geometry of the device. The effect of the device 
dimension is clearly seen in a larger value of Vg^ in device-2 shown in FigJ2]D where the 
doubly clamped beam is 120 nm thick, as against 100 nm for device- 1 in Fig|2^ (the lengths 
differ by 200 nm). 

Next, we consider another complementary aspect of the data - the amplitude of the 
mixing current. FigJSJi shows /^^^.(Aa;) (log-scale) as a function of V^'^*". We see that the 
amplitude of the mixing current for the negative values of are significantly larger than 
those for positive V^^'" for the same mechanical mode; this is also seen in FigJTj. We now 
try to understand this asymmetry as our InAs NW are n-type semiconductors [26| as seen in 
Figllt. To understand this asymmetry we have carried out detailed fits of the experimental 
data for the amplitude of mixing current as a function of frequency using Eqn{T] The 
amplitude z{ijj) of oscillation at frequency co is given by 

^^^^ ^ <- cos(A0 + arctan(^)) ^ 

Q y(i-(^)T + (^)^' 

where is the amplitude at resonant frequency loq, Q is the quality factor and is 
the relative phase difference between the u and u + Au signals that depends on the device 
parameters like the contact resistance. Fitting from EqnsH] and [2] allow us to extract the 
variation of Q as a function of V^'^ as shown in Figj2^ ([ss!]). We have also estimated the 
amplitude of oscillation using Eqns{T] and [2] by examining the ratio /"""^(/^^j = c v (J) — ■ 
A plot of the calculated amplitude (-2^^°) is seen in Figj2]F. We see that as the \V^'-'\ is 
increased Q and zl^° are observed to decrease. We also observe that there are noticeable 



"amp 

differences in Q and for positive and negative values of V^'-' . The values of Q and z^^" 
are larger for positive . One of the possible mechanisms that can explain this behavior 
in the amplitude is that with increasing V^^*" one increases the density of electrons in the 
NW leading to reduction in the screening length. This implies that the simple geometrical 
capacitance is inaccurate particularly since the screening length can be comparable to NW 
diameter at low densities (at negative V^'" in our case). In our device geometry, using 
Thomas- Fermi approximation [22], screening length is around 20-40 nm (diameter of our 
devices are 100 nm) and the distance between the nanowires and gate oxide 200nm. So, 
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the screening length is a significant fraction of the diameter and the suspension distance - 
this plays a critical role in observing the effect of density gradients within the cross-section 
of the nanowire. In case of single walled carbon nanotube, diameter is 1-2 nm and height 



of suspension is typically lOOnm or more 



13l-ll6l|. Additionally the screening length in 

so due to an 



carbon nanotubes is typically several multiples of the nanotube diameter 
increased ratio of suspension distance to diameter and the large screening length compared 
to the diameter it is very difficult to observe the physics we discuss for the case of nanowires 
in carbon nanotubes. We would like to point out that this is not a peculiarity of the InAs 
nanowires and should be seen in other semiconducting nanowire devices as well with similar 

dimensions. Additionally, if one considers the realistic case with non-uniform density of 

I 

carriers in the NW due to the device geometry |37l| the NW will have a gradient of dielectric 
constant jssi]. A gradient of dielectric constant 



alongwith a change in capacitance as a 



function of the density changes the capacitive coupling of the NW to the gate. This results 
in differing amplitudes for two different electron densities. Our device geometry with NW 
diameter comparable to the gap accentuates this effect. 

In order to better understand the effect of the gate voltage in tuning the spatial charge 
density across the cross-section of the nanowire we selfconsistently solve three dimensional 
Poisson's equation using finite element method for our device geometry. We use approach 
followed by Khanal et al. [sT^ by solving V ■ erfV$(x, y, z) = p{x, y, z), throughout the space 
of the nanowire and its dielectric environment (here is dielectric constant, $ is the local 
electrostatic potential in the system due to applied gate voltage, and p is space charge density 
inside the NW). The geometry consists of a 100 nm diameter and 1.5 fim long InAs nanowire 
clamped by metallic electrodes. The wire is suspended 100 nm above a 300 nm thick silicon- 
oxide dielectric on the gate electrode. Inside the nanowire the p{x, y, z) = e(A^d— 
where e is the charge of an electron, A''^ is the density of the n-type dopants (~ lO^^cm"^), 
n($) and are the densities of electrons and holes. The unintentional source of n-type 
dopants in our growth is Si and C, from the metal organic precursors and are assumed to 
be uniformly distributed throughout the nanowire. A selfconsistent calculation gives us the 
distribution of potential throughout the space and space charge density in the nanowire. 
FiglS^ and 3b show the colourscale distribution of potential when the V^^*" = —5 V. In 
order to model the consequences of modifying gate voltage we calculated the distribution of 
space charge density through the nanowire for V^'^ = —25, —5 and 5V. Figj3b, 3d, 3e, and 3f 
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show the result of such a calculation for varying V^'" for space charge density in the center 
of the NW. For positive voltages the Fermi energy is very close to the conduction band of 
the InAs and as a result the charge density is very uniform while behaving as a metal like 
system. This simple modeling supports our arguments that a gradient of electron density 
can modify the capacitive coupling and the resulting amplitude. Further analysis is needed 
to solve self consistent solutions to the Poisson's equation where the dielectric constant is 
itself a function of the density and the quantitative variation of Q and amplitude with sign 
of V^'^ . Our measurements suggest a way to tune the efficiency of actuation by tuning the 
density of carriers. 

We next consider three other features of our data - first, the presence of other mechanical 
modes near the fundamental mode; second, the mixing of modes as a function of the V^^*^, 
and third, the non-linear properties of NW oscillators driven to large amplitudes. FigJl^ 
and 4b show the plot of /^^^.(Aa;) as a function of V^^*" and u for device-3 and device-2. It is 
well known that for a doubly clamped beam with no tension, /„ = Cni2\ — with Cq = 1.78, 
Ci = 4.90 and C2 = 9.63 for the transverse modes. It is clear that if the fundamental mode 
(/o) is described by a mode with zero nodes and moving in a plane perpendicular to the 
substrate, the other observed modes, in the frequency range near /q (Figj4^, 4b), cannot 
be defined by /i and /2. We have observed the anticipated /i and /2 modes at higher 
frequencies. The other modes beside the fundamental in FigJ3^, 4b, are explained due to 
geometrical asymmetry along the diameter in the NW (FigJH:;, details of this calculation 
provided in Supporting Information). These are the modes involving motion in a plane 
that is not perpendicular to the substrate. This would explain the less steep slope of the 
dispersion function of V^'^ . 

FigJl^ and 4b also show the mixing of the modes as a function of V^^*" . The mode mixing 
can be seen clearly in FigJUi which shows a close-up of the data in FigJDD. Displacement 
along the transverse direction y (perpendicular to z) will weakly affect the capacitance 
because of any slight asymmetry in the physical structure of the NW. The coupling coefficient 
|^|^(V^^'")^ appearing in the potential energy gives rise to mode-mixing (see Supporting 
Information). For the device shown in FigJDo the minimum frequency gap in the region of 
level repulsion is 0.25 ± 0.02 MHz. The asymmetry in the amplitude of the modes away 
from the region of mixing can also be understood within this model. 

We next consider the non-linear response of these NEMS oscillators. Due to electrostatic 
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actuation, the potential energy of the oscillator is asymmetric about the equilibrium position 
and has the form V{z) = ^K^ffz"^ + 9z^ + fiz^ + H.O.{z), where K^fj, and fi are coefficients 
depending on V^^*". We have experimentally probed the non-linear and hysteretic response 
of the device. FigJSK, 5b, 5c, and 5d, show the experimentally measured non-linear response 
of this device leading to hysteretic behavior as function of with increasing amplitude of 
driving force. Here, V^'" is swept at a given drive frequency to measure the mixing current 
and several such V^'^ sweeps are collated (the distinct response of sweeping the frequency 
at a fixed Vj^'" - a common mode to study non-linear response - is described in supporting 
information). There are two features that we would like to point out - first, in the region 
at the bottom of the W-shaped dispersion curve, two branches of the same mode merge 
into one broad peak where the oscillator has large amplitude over a wide range of V^'" and 
second, whenever during the V^'" sweep, one crosses the region with a local negative value 
for -^j^uu (here, /o is the resonant frequency at a particular V^'^), shows a curved hysteretic 
j-shaped response (indicated by seen in Fig|5b. In order to understand and explain the 



experimental 



Uy observed hysteretic response as a function of Vg , we have used the Duffing 
equation 40|] for our resonator. The result of such a calculation for amplitude is seen in 
Figj5^, 5f, 5g, and 5h, with increasing amplitude of driving force. To calculate the amplitude 
we have only used the observed dispersion relation (W shape) as input for Duffing equation. 
There is a qualitative agreement between the experimentally measured data shown in FigJS^- 
d and results of calculation using the Duffing equation shown in Figj5fe-h. We find that for 
every increase in excitation amplitude by 100 mV corresponds to an increase of a factor of 
2.5 in the anharmonic component of the Duffing equation (from observing the calculated 
data). Additional aspect of the nonlinearity of the oscillator is also seen in the evolution of 
dispersion near V^'" = as the actuation amplitude is gradually increased from 100 mV to 
400 mV in FigJS^-d. The negative dispersion is due to the softening of electrostatic force 
and with larger amplitude of oscillation the effective spring constant changes, as the wire 
samples a region with varying electric field; this difference is clearly seen in the dispersion 
near V^'" = for the data shown in Figl5] a and d. 

Figj5}, shows the line plot (along the dashed line in FigJSf) of calculated amplitude 
for different sweep direction of V^'" and the resulting hysteresis. The observed hysteretic 
response, as a function of V^'" , is quite different from the hysteretic response as a function 
of drive frequency ( discussed in supporting information ). This non-linear response of our 
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devices with Vg^'" may be utilized for charge detection 2J] as near the onset of non-hnearity 
the change in amphtude as a function of Vg^'" is very large. 

In summary, we have studied the electromechanical properties of doubly clamped InAs 
NW resonators. Their size and tunable electron density allow us to map behavior that has 
not been manifested in a single device. We have observed and quantitatively explained the 
competition between the softening of stiffness of the restoring force of the resonator due to 
the variation of the electrostatic force under variable gate voltage. At larger voltages, the 
stretching of the nanowire leads to increased stiffness resulting in a non-monotonic dispersion 
of the fundamental mode with V^'" . The screening of electric fields due to the variation in 
the density of the electrons in our suspended FET devices modifies the amplitude because 
the variation of the screening length spans the cross sectional dimension of our nanowire. 
Further, the non-linear properties of our device can be understood qualitatively using the 
Duffing equation that explains the hysteretic response of the amplitude as a function of 
gate voltage. Thus in a single device, we demonstrate, separate and account for three 
diverse behaviors. Our measurements indicate that measuring electromechanical response 



influenced by charge screening could lead, in the 



uture, to new ways to probe spin physics 
by exploiting spin-dependent charge screening [4l|]. Probin g th e physics by tuning electron 
density in NEMS devices may help probe the role of defects 42|] and electron hopping as one 
moves from insulating to conducting regimes. Control over the non-linear dynamics may be 
achievable by controlling mode mixing and DC gate voltage. 

This work was supported by Government of India. J. Parpia was supported by DMR- 
0457533. 
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FIG. 1: (color online) a) Tilted angle SEM image with the circuit used for actuation and detection 
of resonance for an InAs NW resonator. The diameter of the wire is 100 nm and the length of NW 
is 2.9 fj,m. The scale bar indicates a length of 200 nm. b) The mixing current (Aoj/lir = 17 KHz) 
as a function of frequency for two values of V^'^' (mixing current shown for +ve V^'^ , is 5 times of 
its original value), c) Variation of the conductance as a function of DC gate voltage, d) The plot 
of sound velocity Vs = calculated using the measured frequency {uj/2tt) of the fundamental 
mode of the InAs NW resonators and geometrical values. The dashed line indicates the speed of 



sound obtained for Vg = \/E/p using bulk values for E and p 
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FIG. 2: (color online) a) and b) Color logscale plots of mixing current as a function of V^'-' and 
uj/2tt for two devices (for device-1, Figl2^, diameter d=100 nm, length 1=2.9 /im, for device-2, 
FigEb, d=120 nm, 1= 3.1 fim). c) Calculated dispersion as a function of ^^'^ ■ d) Lineplot at 29 
MHz for device-1 (dashed line in Fig|2^). e) The plot of the Q as a function of V^^*^ for device-1. 
f) The plot of the "amplitude" ^^mp ^ ^ function of V^'^ for device-1. (The red and blue traces 
in Figj2^, and in Figj2f, show the data for negative and positive gate voltages). 
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FIG. 3: (color online) a) and b) Shows the result of a FEM based self-consistent solution of Poisson's 
equation giving the potential around a doubly clamped suspended nanowire device 100 nm in 
diameter and 1.5 /um long. The separation of the NW and 300 nm thick Si02 is 100 nm. The 
back plane of the nanowire is the gate held at -5V and the two terminals of the wire are grounded. 
The maximum of the colorscale bar (red)is OV and the minimum (blue) is -5 V. c) The space 
distribution in the cross section of the nanowire for = -25V shows the gradient. The log- 

colorscale below varies from 10 C/m^ (blue) to 10^ C/m'^ (red), d) Plot of space charge density in 
the vertical direction at V^^*^ = -25V through the middle of the wire. The asymmetry along the 
vertical direction due the gate below the wire is clearly seen, e) The space distribution in the cross 
section of the nanowire for V^'^ = -5V shows the gradient. The log-colorscale below varies from 
10^ C/m^ (blue) to 10^ C/m^ (red), f) The space distribution in the cross section of the nanowire 
for Vg= 5V shows uniform space charge distribution as the electron density in the nanowire is 
increased when the nanowire field effect transistor is turned on. The log-colorscale below varies 
from 10^ C/m^ (blue) to 10^ C/m^ (red). 
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FIG. 4: (color online) a) and b) Shows mode mixing for device-3 and device-1 respectively (for 
device-3, d=103 nm and 1=3.1 ^m). c) Calculated dispersion for mode mixing using asymmetry 
in the wire, d) Zoomed-in view for device-1 (FigUja) which shows the mixing of the modes as a 
function of ■ Minimum separation between the modes, indicated by the double-headed arrow, 
is 0.25 ± 0.02 MHz. 



15 




-40 -20 20 40 -40 -20 20 40 
Gate Voltage (V) 




Gate Voltage (V) 

FIG. 5: (color online) a), b), c), and d) Color logscale plots of mixing current from device-1 as 
a function of V^^*^ and uj in non-linear regime with increasing amplitude of driving force starting 
from 100 mV to 400 mV respectively, e), f), g), and h) Colorscale plots of calculated amplitude of 
a resonator in the non-linear regime using Duffing equation with increasing amplitude of driving 
force from Figj5^ to Figj5)i. i) Shows the amplitude for the two directions of V^'-^ which indicates 
hysteresis at constant frequency 31 MHz (red and blue arrows show corresponding gate sweep 
direction). 
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